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Abstract 

Often noisy point clouds are given as an approximation of a particular compact set of interest. 
A finite point cloud is a compact set. This paper proves a reconstruction theorem which gives a 
sufficient condition, as a bound on the Hausdorff distance between two compact sets, for when 
certain offsets of these two sets are homotopic in terms of the absence of //-critical points in an 
annular region. Since an offset of a set deformation retracts to the set itself provided that there 
are no critical points of the distance function nearby, we can use this theorem to show when 
the offset of a point cloud is homotopy equivalent to the set it is sampled from. The ambient 
space can be any Riemannian manifold but we focus on ambient manifolds which have nowhere 
negative curvature. In the process, we prove stability theorems for /i-critical points when the 
ambient space is a manifold. 

1 Introduction 

In modern science and engineering a common problem is understanding some shape from a point 
cloud sampled from that shape. This point cloud should be thought of as some finite number of 
(potentially) noisy samples. Topology and geometry are considered very natural tools in such data 
analysis (see e.g. [H] and [1]). One reason is because topological invariants are often more stable 
under noise. We will want to understand the homotopy type of a set - two objects are homotopy 
equivalent if there is a way to continuously deform one object into another. 

Often we wish to know the extent to which we can, and how to, reconstruct shapes from noisy 
point clouds. Naturally the more restrictions on the space the easier it is to reconstruct. The first 
area of focus was the study of surfaces in M 3 , motivated by problems such as medical imaging, 
visualization and reverse engineering of physical objects. Algorithms with theoretical guaranties 
exist for smooth closed surfaces with sufficient dense samples. The sampling conditions are based 
on the concept of e-sampling. A point cloud of A is an e-sample if for every point p £ A there is 
some sample point at distance at most elfs(p) away, where lfs(p) (called the local feature size at p) 
is the distance from p to the medial axis of A. The Cocone Algorithm produces a homeomorphic set 
from any 0.06-sampling[l] of a smooth closed surface. This process has been extended to smooth 
surfaces with boundaries [9]. However given an arbitrary compact set, the best we can hope for is 
to find some nearby set that is homotopy equivalent. For instance from a point cloud we can not 
tell apart the original set and the same set with a slight thickening in places. 

One of the simplest methods of reconstruction is to use the offset of a sampling. Given a set K, 
the r-offset of K, denoted K r , is defined to be {x £ M. : d{x,K) < r}. This is topologically the 
same as taking the a-shape of data points |12| or taking the Cech complex [6]. This leads to the 
problem of finding theoretical guarantees as to when an offset of a sampling has the same topology 
(i.e. homotopy type) as the underlying set. In other words we want to find conditions on a point 
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cloud S of a compact set A so that S r is homotopy equivalent to A. We will in fact find sufficient 
conditions for S r to deformation retract to A. Clearly this will only work if the point cloud is 
sufficiently close. Usually "sufficiently close" is interpreted as a bound on the Hausdorff distance 
between A and S (the Hausdorff distance between A and S, denoted dn(A, S), is the smallest r > 
such that S C A r and A C S r ). Much earlier theory assumes that this Hausdorff distance is less 
than some measure of geometrical or topological feature size of the shape and show the output is 
correct. We now survey some of this development. 

The medial axis of a compact set A is the set of points p in the ambient space for which there 
is more than one point in A which is closest to p. The reach of A is the minimum distance between 
points in A and points in the medial axis of A. It can be thought of as the minimum local feature size. 
Sampling conditions based on reach include those found in [20] which consider smooth manifolds in 
W 1 . Smooth submanifolds have positive reach but a wedge, for instance, does not. To deal with a 
larger class of sets Chazal, Cohen-Steinerm and Lieutier in [7] introduced the notion of //-reach. A 
point is //-critical when the norm of the gradient of the distance function at that point is less than 
or equal to // (a geometric description is expanded on later - see section 3). In particular 0-critical 
points are critical points of the distance function and every point on the medial axis is a //-critical 
point for some // < 1. The //-reach of a set is the supremum of r > such that the r offset such not 
contain any //-critical points. In [7j and [2] sampling conditions are given in terms of the //-reach. 

Instead of making our bounds in terms of //-reach, we will only require the absence of //-critical 
points in an annular region of A along with a bound on the weak feature size. The weak feature 
size is the infimum of the positive critical values of the distance function from A. This has several 
advantages. Firstly, it can mean a significant improvement on the bounds such as in the case of 
"hairy" objects. Secondly, it can be applied to a larger class of compact sets. Instead of positive 
//-reach, all that is necessary for our reconstruction process to work is for A to have positive weak 
feature size. Compact subsets with a cusp or including two tangential spheres are examples of 
compact sets with zero //-reach for all // > but these still have positive weak feature size. Every 
semi- algebraic set has positive weak feature size |13j . This follows from the fact that the distance 
function from a semi-algebraic set has only finitely many critical values. 



One limitation to any of these reconstruction theorems is the requirement of knowing geometric 
properties of the unknown object we are trying to reconstruct. A shift in perspective can overcome 
this limitation by considering the geometric properties of the point cloud, which we do know, and 
can hence prove sufficient conditions for an offset of the original set to deformation retract to an 
offset of the point cloud. We know the point cloud and hence we know the //-critical values of its 
distance function. Theoretical guarantees are given in [7J for when suitable homotopies exist by 
considering the /i-reach of an offset of the point cloud. Unfortunately there is usually a significant 
number of small critical values of the distance function to the point cloud. This means the starting 
offset beyond which //-reach is considered is significant. Our approach which only considers the 
existence of //-critical points in a annular region thus gains a significant advantage. 

We note that previous reconstruction theories have been restricted to the case where the ambient 
space is Euclidean. Another contribution of this paper is to allow the ambient space to be any 
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manifold whose curvature is bounded from below, thus answering an open question asked in [7J. 
Although we focus on the important case of non-negatively curved manifolds we explore a paradigm 
of reconstruction which can be applied to manifolds with curvature bounded from below by some n < 
with analogous, albeit messier, results. Examples of manifolds with nowhere negative curvature 
are Stiefel and Grassmannian manifolds. These examples are important because there are many 
applications where data naturally lies on these manifolds such as in dynamic textures |10j . face 
recognition [5], gait recognition [3] and affine shape analysis and image analysis [21]. 

Even when restricted to the case where we use //-reach in Euclidean space we still improve on 
the previous results whenever // < 0.945. Our main theorem expressed as a sampling condition is 
as follows. 

Theorem. Let fi G (0, 1), r > 0. Let A be a compact subset of a smooth manifold Ai with nowhere 
negative curvature such that the injectivity radius of every point in A4 is greater than r and A r is 
homotopic to A. Let S be a (potentially noisy) finite point cloud of A (i.e. a finite set of points). 
Suppose that either 

(i) there are no fi-critical points of the distance function from A in {x G M : d(x, A) G [a, b]} 



(ii) there are no fi-critical points of the distance function from S in {x G M. : d{x, S) G [a, b]}. 



Note that a sufficient condition for A a to defamation retract to A is that the weak feature size 
of A is at least a [E] . 

Our approach is to show the existence of a Lipschitz unit vector field whose integral flow induces 
a deformation retraction from S r to A r ^s or from A r to S r -s- For the flow of a vector field to have 
a certain desirable property, we consider sufficient local conditions on potential unit tangent vectors 
at each point. To this end we define cone fields; cone fields are generalizations of not necessarily 
continuous unit vector fields where we attach a closed ball in the unit tangent sphere, a "cone", 
at each point in the manifold. To prove the existence of useful cone fields we need to generalize 
previous stability results about //-critical points. A stability result of //-critical points when the 
ambient space is Euclidean is proved in [7J. We prove a generalization of this result for when the 
curvature of the ambient space is bounded from below. It is worth observing that although //-reach 
is not stable under Hausdorff distanc^] we do have some stability of an absence of //-critical points 
within of annular regions. A quantitative analysis of this phenomenon leads to the theorem. 

2 Notation and cone fields 

We first need to establish some notation. Throughout (A4 , g) is a smooth n-dimensional manifold 
without boundary and K C M. is a compact subset. We will denote the closed ball of radius r and 
centre x by B(x,r). Let dx denote the distance function from K. We will use for convenience: 
K$ := (f^fO,^ and K\ a ,b] := d~^~[a,b\. The Hausdorff distance between two compact sets K and L 
is denoted dn(K, L) and is defined by 



1 In particular, for any compact set K and any bound on Hausdorff distance S > there is a compact set L with 
zero /i-reach such that cLh {K, L) < 5 



or 



Then S r is homotopic to A whenever dn{S, A) < min < r — a 
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dn(K, L) := maxjsup di(x), sup dxiy)}- 
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Alternatively it is the smallest r > such that K C L r and L C K r . 

Let K° denote the interior of K. Given a metric space X, let KL{X) denote the metric space of 
compact subsets of X with the metric determined by Hausdorff distance. 

We want to recall some differential geometry. A useful reference is [T7] . The unit tangent bundle 
of manifold a (M,g), denoted by UTAi, is the unit sphere bundle for the tangent bundle TAi. It 
is a fiber bundle over At whose fiber at each point is the unit sphere in the tangent plane; 

UTM :=]J{ve T X (M) : g x (v,v) = 1} , 

where T X M denotes the tangent space to M at x. Thus, elements of UTM are pairs (x, v), where 
x is some point of the manifold and v is some tangent direction (of unit length) to the manifold at 
x. 

The exponential map at x is a map from the tangent space T X M to M. For any v G T X M, a 
tangent vector to M at x, there is a unique geodesic j v satisfying 7«(0) = x with initial tangent 
vector 7^,(0) = v. The exponential map at x is defined by exp x (v) = 7 V (1). The injectivity radius 
at a point x is the radiusof the largest ball on which the exponential map at x is a diffeomorphism. 
Normal coordinates at a point x are a local coordinate system in a neighborhood of x obtained by 
applying the exponential map to the tangent space at x. 

Consider the (n — l)-dimensional unit sphere, S n , lying inside R n with a metric induced from 
this embedding. Denote by C(w,/3) the closed ball in S" 1-1 centered at v with radius ft. We can 
view C(w,f3) as the intersection of S" 1-1 with a particular infinite cone; 

C(w, (3) = S n - X D{ve R n \{0} : Z(v, w) < (3}. 

Here all the angles are in radians. We say C(w,(3) is acute if f3 is acute. 

Denote by F the fibre bundle over M where each fibre over x G M is the space of non-empty 
closed balls in the unit tangent sphere at x. F has a natural metric induced from the Hausdorff 
metric on compact subsets of UTAi. 

A cone field over a subset U C M is a section of F restricted to U. A cone field is continuous if 
it is continuous as a section. We can consider a cone field as a section from U C Ai into JC(UTM), 
the space of compact subsets of UTAi with Hausdorff metric, and it is continuous if this map is 
continuous. As a set, we can write cone fields as 

{(x,C(w x ,P x )) : w x G UT X M,P X 6 [0,vr]}. 

Given two cone fields V and W we will write W > V to mean that for any point x the cone at x in 
V is contained in the cone at x in W. 

Let K be a compact subset of W 1 not containing the origin. We say C(w,f3) is a minimal 
spanning cone for K if K := {A : y £ K} C C(w,(3) and K C C(w',(3') implies /3 < f3' . 

Lemma 2.1. Let K C W 1 be a non-empty compact subset set such that is not in the convex hull 
of K . Then there exists a unique minimal spanning cone of K which is acute. Furthermore the 
map from compact subsets of M rf which don't have the origin in their convex hull to their minimal 
spanning cone is locally Lipschitz. If we restrict our attention to compact subsets whose distance to 
the origin is bounded from below by r we further have that the map sending each compact subset to 
its minimal spanning cone is Lipschitz with Lipschitz constant bounded from above by 1/r. 

Proof. From the compactness of K we know that, given any v G S"™ -1 , there is some angle j3 v such 
that K C C(v, (3 V ) and if K C C(v, (3' v ) then /3 V < /3' v . Observe that for many choices of v here the 
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corresponding (3 V is not acute. Set j3 := inf{/3„ : v G S*™" 1 }. The map v >->■ /3 V is clearly continuous 
and hence the infimum is achieved. Any C(t> , f3 v ) where /3 V = /3 is by definition a minimal spanning 
cone for i^T. 

Since is not in the convex hull of K, there is some set of coordinates such that K lies in the 
upper half plane. This implies that is some w such that K C C(w, 7r/2)°. This implies that the 
corresponding angle (3 W is acute. This means that any minimal spanning cone of K is acute. 

Thus all it remains to show is that there is only one w G S 11 ^ 1 such that f3 w = (3. Suppose that 
there are two distinct such vectors w\ and W2- Now K C C(wi,f3) and K C C(w2,f3) and hence 
K C C(wi,(3) n C(w2,(3). The angle between w\ and W2 must be less that 2(3 < it as otherwise 
C(wi,P) n C(u>2, (3) would be empty. Let 9 = Z(ti>i, u^)- Set := (wi + 102) /\/2 + 2cos# which is 
a unit vector by construction. 

Let x £ K and u the unit vector in the direction of x from 0. We have Z(v,w\) < (3 and 
/-(y , W2) < /3 but since w\, W2 and t> are all unit vectors we have (v, wi) > cos (3 and {v, W2) > cos (3. 
From the choice of w we have (v,u>) > 2 cos (3/^2 + 2 cos # which implies that Z(-u,u>) =: a x is 
strictly less than (3. 

Since K is compact we have a := supja^ : x G i^T} is attained and a < /3. However this implies 
that K is contained in a cone with an angle smaller than (3 which is a contradiction. Thus the 
choice of centre vector w is unique. 

Let K be a compact subset of M n whose convex hull does not contain the origin with a corre- 
sponding spanning cone C(K). There is some 5 > such that K28 also does not contain the origin 
in its convex hull. Let L be a compact subset of W 1 such that dn(K,L) < e < 5. Let ir be the 
canonical projection map from M n \0 onto S^ 1 . We have cIh{tt(K),tt(L)) < dn{K, L)/r < e/r. 
where r = min{||x|| : x G -ftT^} > 0. 

Since LcKgwe have that the convex hull of L does not contain the origin and hence there exists 
a minimal spanning cone C(L) furthermore L e C K25 and hence L e also has a minimal spanning 
cone C(L e ). 

Now vr(L) C ir(K e ) C C(ir(K e )) C C(vr(K)) e/r implying C(L) = C(vr(L)) C C(K) e/r and 
similarly we know C(ZT) C C(L) e/r . Together they imply that d H (C(K), C{L)) < e/r. □ 

We can consider vectors inside the cone at a point x. We say a vector field X := {(x,v x );x G 
U, v x G T X M} is subordinate to the cone field W = {(x, c x )} if v x always lies in c x . We will call a 
vector field strictly subordinate if the vector at each point lies in the interior of the cone. 

Given a compact subset K C Ai, a point x G Ai\K and some R > 0, consider the compact 
subset exp" 1 i^ni?(0, i?) of T x Ai. The straight line segments from to points in exp" 1 K (1B(0, R) 
are sent under exp x to geodesies of length at most R from x to points in K. If is not in the 
convex hull of exp" 1 K n B(0, R) then by the previous lemma there exists a unique acute minimal 
spanning cone for exp" 1 ^!"! B(0,R) from 0. Considered as a subset of UT x Ai, we call this the 
minimal spanning cone for K from x of length R. If, for some subset U G Ai, is not in the convex 
hull of exp" 1 K n B(0, R) for any x G J7 then we can define over [/ the minimal cone field for K 
with length by assigning to each x the corresponding minimal spanning cone. 

In general, even if a minimal spanning cone field exists over U G Ai it may not be continuous. 
However we have the following sufficient condition for the minimal cone field to be locally Lipschitz. 

Proposition 2.2. Let R > 0, and K, L compact subsets of Ai such that the injectivity radius of 
every point in K is greater than R, K C B(x, R) for all x G L, and K n L = 0. If there exists an 
acute minimal cone field over L for K of length R then that minimal cone field is Lipschitz. 

Proof. Since K and L are compact and K (1 L = there is some e > such that L t n K = and 
the injectivity radius for every point in K is at least R + e 
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For each fixed y E K we can construct the map 

where Vx is the tangent vector to the unique shortest geodesic from x to y with length d(x, y). This 
is well-defined because of our injectivity radius assumption and because K n L = 0. The map y is 
smooth. Since L is compact and L <Z L° we can conclude that <^|l is Lipschitz with some Lipschitz 
constant a y . 

Set a := sup{a^ : y E if} which is finite as K is compact and the map y \- > a y is continuous. 
This means every (f> y is Lipschitz with Lipschitz constant bounded from above by a. 
Now consider the map 

: L -> K(TM) 
x ^ (J vg. 

Observe that dni^yix), 4>{x')) < dn{<i>y{x), <f> y (x')) and dii{(l>(%),(f>y{x')) < dif((f>y{x),<i>y(x')) for all 
y £ K. This implies that 

c? fl -((/>(x),0(x / )) < sup{dH(0j / (^),0y(2; / )) : y e if} < ad(x,x'). 

and hence is also Lipschitz. 

From y) > e for all x 6 L and y G if we know that (2(0, Uj/eK > e - R- ecan that in 
Lemma |2.1| we showed that the map that sends compact subsets whose distance to the origin is 
bounded from below by e and do not have in there convex hull to their minimal spanning cone 
is Lipschitz. Since the composition of Lipschitz maps is Lipschitz, we conclude that the minimal 
spanning cone field is Lipschitz. □ 

We can define the complementary cone for all acute cones. The complementary cone of C(w, (3) 
is C(w,ir/2 — f3). Given a cone field where the cone at each point is acute we can construct the 
complementary cone field pointwise. From the triangle inequality on the unit tangent sphere we 
obtain the following useful lemma. 

Lemma 2.3. Let v be a vector in some acute cone C and v a vector strictly inside the comple- 
mentary cone. Then Z(v,v') < n/2. 

3 Stability of ^-critical points 

We want to study the gradient vector fields for distance functions from compact subsets of a general 
manifold (M,g). This can be thought of as the obvious generalization of the gradient vector fields 
for distance functions from compact subsets of Euclidean space (as studied in [7] and [TH]). 

We call 7 a segment from x ^ K to K is if 7 is a distance achieving path from x to K. If A4 is 
Euclidean then these segments are straight lines. Observe that on an arbitrary manifold there can 
be more than one segment connecting x to the same y £ K. The point a; is a called a critical point 
of the distance function from K if, for all non-zero v E T X A4 , there exists a segment 7 from x to K 
such that Z(7'(0),i>) < tt/2. Equivalently, if T(x) := {y E K : djc(x) = d(x,y)}, then x is a critical 
point if and only if lies in the convex hull of exp" 1 T(x) H -6(0, dx(x)). 
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We need to construct the gradient vector field so that it vanishes at critical points of the distance 
function. For all non-critical points we can consider the minimal spanning cone C(w x ,{3 x ) for T(x) 
from x of length dx{x). We set 

Vk(x) := -cos((3 x )w x 

whenever x is not critical. Observe that Vx a (i) = V^(x) whenever dx(x) > a. For [i G M, we 
called x ^-critical if [|Vjf(ac)[| < \i. A point is O-critical exactly when it is a critical point for the 
distance function. 

It is easy to verify that these definitions agree with those given in [7J when Ai is Euclidean. 

We will want to prove a generalization of the stability result in [7] where the ambient space is 
a manifold with curvature bounded from below in a suitable neighborhood of the compact subset 
under study. By appropriate scaling it is sufficient to consider the cases where the curvature is 
bounded from below by 1, or —1. 

Lemma 3.1. Let K C M be a compact subset of a manifold (Ai,g) for which the curvature on 
is bounded from below by k. Let x G K a \K be a ^-critical point of dx- Then for any y G K^a 
we have 

cos dxiy) > cos d(y, x) cos dx{x) — sind(y, x) sindx {x)n if k = 1 

dxiy) 2 <dR(x) 2 + d(x, y) 2 + 2d(y, x)d K (x)n if K = 

coshdx(y) < cosh d(y, x) cosh dx(x) + sinhd(y, x) sinh if n = — 1. 

Proof. Let 8 G (0,7r/2] such that cos 8 = fi. Let T(x) = {z £ K : dx(z) = d(z,x)} and set 

f^):={z/\\z\\:zGe W -\r(x))}. 

Fix y £ K2a and choose v G T x Ai such that exp x (d(x,y)v) = y. 

We want to show that there is some z £ K and length achieving geodesies j x and ^y x such that 
Z(y,x,z) < ir — 8 where Z(y,x,z) is the angle between 7^ and 7*. Suppose not. This means that 
no point of T(x) lies in C(v , ir — 8). Geometrically this means that T(x) must lie in C(—v, 8)° which 
is the complement of C(v,tt — 8) in the sphere (recall that C(—v, 8)° is the interior of C(—v, 8)). 

However this implies that the minimal spanning cone for T(x) from x of length dx(x) lies strictly 
inside C(—v,8) and hence ||Vif(x)|| > cos 8 = fi. This contradicts the assumption that x is a [i- 
critical point (i.e. ||Vif(x)|| < [i). Thus by contradiction, there is some point z G K and length 
achieving geodesies 7^ and 7^ such that Z(y, x, z) < tt — 8. 

Let ^ x , y ,z be the geodesic triangle with 7^ and 7! such that Z(y, x, z) < ir — 8 which must exist. 
Let Ax,j/,i be the corresponding triangle in A4(k), the manifold with constant curvature k, where 
the length of the sides are preserved. Toponogov's theorem is a triangle comparison theorem which 
quantifies the assertion that a pair of geodesies emanating from the same point spread apart more 
slowly in a region of high curvature than they would in a region of low curvature. The details of this 
theorem and its proof can be found in [Sj. By taking the contrapositive of Toponogov's theorem we 
know Z(y, x, z) < Z(y, x,z) < ir — 8 and hence cos Z(y, x, z) > — [i. 

We finally substitute dxiy) < d(y, z), dx{x) = d(x, z) and cos Z(y, x, z) > —\x into the spherical, 
Euclidean and hyperbolic cosine rules respectively to obtain the desired inequalities. □ 

Our stability result will arise from comparing two opposing inequalities - one from the previous 
lemma alongside one coming from the following lemma which is Lemma 4.1 in |19j . It is easy to 
check the definition of ||Vif (x)|| coincides with ||V X /|| when / = dx and X is a manifoldj^] 

2 11^/11 is the nonnegative number max{0, limsupy^^. ^^~{Jj x ' > }■ That this is ||Vif(ic)|| follows from our geometric 
construction of Vk{x) and from the cosine rule. 
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Lemma 3.2 (Lemma 4.1 in |19|). Let X be a metric space. Suppose f : X — >■ M is a locally 
Lipschitz map, x G X, and /(x) = 0. For /i, r > ; assume that the ball B(x,r)) is complete and 
that || V z /|| > /i /or each z with d(z, x) < r and f(z) > 0. Then for each < C < /J, there is a point 
z G X with d(z,x) < r and f(z) = Cr. 

The following proposition is a generalization of critical point stability theorem in |7] where the 
ambient space is can now any manifold with non-negative curvature. 

Proposition 3.3. Let K, L be compact subsets of M. with dn(K, L) < 6. Let x be a \i-critical point 
ofd K . If 

C > n + 



d K {x) 

and K dK ^ +i s/(c-/i) has nowhere negative curvature, then there exists a C-critical point y of d^ 
with dL{y) > di(x) and y G B(x, 46/ (C — n)) 

Proof. We want to show that there is some y such that ||Vl(x)|| < C and di(y) > di,{x) and 
d(x,y) < 46/{C — (J,). Suppose not. Then there is some p, > C such that ||Vl(x)|| > jl whenever 
dh(y) > d L {x) and d{x, y) < 46/ (C - //). 

If C > n + 2^6/d K {x) then d K (x) - 45/{C - /i) 2 > and hence we can construct 

K '■= K d K {x)-A5/{C-ii) 2 and L := L dK ( x - ) _ 4S /( C _^2. 

By construction dn{K' , V) < dn{K, L) < 6 and d,K'(x) = 46/{C — [i) 2 . Using / = dy — dy{x) and 
r = 46/ (C — //) in Lemma 3.2 we know that there exists a point y G B(x,46/(C — fi)) such that 
f(y) = C46/{C - (j,) which means d L >(y) = d L >(x) + C46/(C - fi). Using d H (K',L') < 6 we can 
show that 

d K '{y) > d v (y) -6 = d L ,{x) + C46/(C - fj) - 6 > d K ,(x) + C46/(C - n) - 26. 

Since dx'(x) = 46/{C — [i) 2 we conclude that 

A5 + CA5(C-y)-25{C- l if 
d K .{y)> jc^) 2 • (1) 



At the same time, Lemma 3.1 implies that dx'iy) 2 < dK'{x) 2 + d(x,y) 2 + 2d(y,x)dK'(x)fM and 
hence 

2 16^ 2 + 16^(C-/i) 2 + 32 ( 5V(C-/i) 

^ (c^ 5 ( } 

By combining ([I]) and ^ and performing some algebraic manipulation we obtain 

1 + (C - /i) 2 + 2(C - > (1 + C(C - /i) - (C* - /i) 2 /2) 2 - (3) 

However algebraic manipulation of ([3]) implies > (C — n) 2 /4 + C[i which is a contradiction. □ 

One of the problems with working with the //-reach is that it is not stable under Hausdorff 
distance. Indeed by the creation of an arbitrarily small cusp we know that for any compact subset 
K, and any 6 > 0, there exists some compact subset L with dn(K, L) < 6 whose ^i-reach is zero for 
all \x > 0. However by only considering /i-critical points in an annular regions we can have stability 
results. 
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Corollary 3.4. Let K, L be compact subsets of a manifold with non-negative sectional curvature 
such that dn{K,L) < 5. Suppose that there are no C-critical points for di in the annular region 

L [a,b]- If 

C>n + 2 



a + 5 

then there are no fx- critical points for dx in the annular region K[ a +s,b-4:5/(C-n)-S]- 



Proof. If x is a /i-critical point with dx(x) 6 [a + 5, b — A5/{C — /x) — 6] then by Proposition 3.3 
there exists some C-critical point y with 

d L (y) 6 [d L (x),d L (x) + 4<5/(C - »)} C [d K (x) - 5, d K {x) + 45/(C - //) + 5} C [a, b] 

which is a contradiction. □ 

Analogous stability results should hold for the cases when k = 1, —1. However, we will only later 
require the case when /i = and so to significantly simplify calculations we restrict to this case. 

Proposition 3.5. Let K,L be compact subsets of Ai with dn(K,L) < 5. Let x be a critical point 
of dx- Suppose that the sectional curvature of K2d K ( x ) ^ s bounded from below by k = —1. 

Then for all C > there exists a C-critical point y of di with d^{y) > di{x) and d(x, y) < A5/C 
whenever 

95 < 2tanh(d K {x))C 2 . 

Proof. Let C £ (0, 1). We want to show that there is some point y such that ||Vz,(y)|| < C and 
d>L(y) > d>L(%) and d(x,y) < 45/C. Suppose not. Then there is some ft > C such that ||Vi(y)|| > jl 
whenever dL{y) > di{x) and d(x,y) < 45/C 

Using f = — dL{x) and r = 45/C in Lemma 3.2 we know that there exists a point y £ 
B(x,45/C) such that f(y) = AS; i.e. d L (y) = d L {x) + A5. From d H (K,L) < 5 we know that 
d K {y) > d K {x) + 25. 

At the same time, Lemma |3.1| implies that cosh cift-(y) < cosh d(y, x) cosh dx{x) and hence 
cosh{d K {x) + CA5/C -25) < cosh(45/C) coshd K (x). 

Using the hyperbolic cosh sum formula and dividing through by cosh dx(x), this can be rewritten 

as 

cosh(45/C) - cosh(2<5) > tanh(d K {x)) sinh(25). 
Our assumption that 95 < 2 tanh(d^(x))C 2 implies that A5/C < 8tanh(dK(x))C /9 < 1. Now 

9 

cosh(i) - cosh(2<5) < cosh(t) - 1 < — r 

16 



whenever t € (0, 1) and sinh(t) > t for all t. This means that we can conclude that 

> tanh(d^(x))2<5 



9 /45 x 2 



16 

and hence that 95 > 2 tanh(<ix(^))C 2 . This contradicts our assumption of 5 implying that there 
must exist a suitably nearby C-critical point. □ 
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4 Reconstruction theorem 



Our reconstruction proof will involve finding sufficient conditions for when there exist minimal cone 
fields. We will use the stability of /i-critical points to show the existence of acute minimal spanning 
cones of As from points in an annular region. 

Lemma 4.1. Let C G (0, 1) and R > 5 > and M be a manifold with nowhere negative curva- 
ture. Let K C M be a compact subset and x E Kr+s- If there are no C-critical points of dp: in 

K [d K (x),d K (x)+2(R-d K (x)+S)/C] and 

4-C 2 
5 < d K {x) - JT ^ J R 

there is an acute minimal spanning cone for Kg from x with length R. 

Furthermore, for a G (0, R + e), this means that if there are no C-critical points of dx in 

K [a,a+2(R+8-a)/C] and 



6 < a 



C 2 



;R 



4 + C 2 ' 

then there is a minimal cone field over K^ a ^ R+ ^ for Kg with length R. 

Proof. Suppose, by way of contradiction, that is in the convex hull of exp" 1 K$ n B(Q, R). 
Set K := (exp" 1 K s ) n B(0, R). Define the map 

R+(d K (x)-5Y 



ip : B(0,R) -»■ dB 0, 



Z H> 



R + (d K {x)-5) z 



and set L to be <p{K). By construction dt{x) = \{R + (dj{ (%) — 5)) and 



d H (exp x L,exp x K) < ^(R- {d K (x) -5)). 



exp x l {K 5 








Construction of L = ip{K) 

Set L = exp x L U {a G K : d(a, x) > R}. Since taking the union with both sets by the same set 
can only decrease the Hausdorff distance, dn{K-,L) < dn(exp x K, exp x L). Also, by construction, 
d L {x) = d v (x) = \(R + (d K {x)-5)). 
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Since is in the convex hull of K (by assumption), there exists z%, . . . , z m G A and oj, . . . a m > 
such that Ti^L^iZi = 0. However this gives 

and hence is in the convex hull of L. Since all the points in L are equidistant from 0, and hence all 
the points in exp L are equidistant to x, we conclude that x is a 0-critical point of d ? . Adding 

exp^ Li 

points further away from x leaves this quality invariant and hence x is a 0-critical point of dx- 
Our condition that 5 < dx{x) — R(A — C 2 )/(4 + C 2 ) can be rewritten as 

\ (R - (d K (x) -8))<~(R+ (d K (x) -S))°^. 

This implies that d H (K,L) < d L (x)C 2 /A, or in other words C > + 2 v / d H (K, L)/d L {x). This 
enables us to apply Proposition |3.3| to say that there exists a C-critical point y of dx with 



d K (y) e [d K (x),d K (x) + 4d H (K,L)/C] = [d K (x),d K (x) + 2(R-d A (x) + 5)/C]. 

This contradicts our assumption about the absence of C-critical points in that annular region. 

For the second part of the proposition we only need to observe that if x G K[ a ,R+S\ then 
[d K (x),d K (x) + 2(R-d K (x)+5)/C]c[a,a + 2(R-a + 5)/C]. □ 



The process can be applied for the case when k = — 1 using Proposition 3.5 instead of Proposition 



3.3 This leads to the following lemma. 



Lemma 4.2. Let C G (0, 1) and R > 5 > and M be a manifold with curvature bounded from 
below by k = —1. Let K C Ai be a compact subset and x G Kr+s- If there are no C-critical points 
of d K in K[ dK ( x ) idK ( x ) +4S/C ] and 

9(R + 5- d K {x)) < 4tanh ( i(H -5 + d K (x)) C 2 



2 

there is an acute minimal spanning cone for K$ from x with length R. 

We observe that so far we have only proved, under certain conditions, the existence of a minimal 
cone field. Note that the minimal cone field is unlikely to be continuous. What we will need is the 
existence of a Lipschitz vector field subordinate to it. The importance of minimal cone fields is that 
vector fields subordinate to them can build deformation retractions. 

Lemma 4.3. Let < 5 < r, and let K, L be compact subsets of a manifold Ai with djj(K, L) < 5. 
Suppose that there exists an acute minimal cone field W over K\ r _g ir+ s] f or K& with parameter r. 
Let W be the complementary cone field to W. If X is a Lipschitz vector field strictly subordinate 
to W' , then X induces a deformation retraction from L r to K r _s- 

Proof. Since X is a Lipschitz vector field it has a unique continuous integral flow (a standard result, 
for example |16j). The idea is to follow this flow from each point in L r until it reaches K r -$. 



Denote the minimal cone field W pointwise by {(x, C(w x , (3 X ))}- From Lemma 2.3 we know that 
any vector sitting strictly inside C(w x , \ — f3 x ) forms an acute angle when paired with angle vector 
in C(w x ,p x ). 

Let X be a vector field strictly subordinate to W' . Let x G L r n K[ r _s^ r+ s] > an d let v be the 
vector at x in X. Now x G B(y, r) for some y G L. Let j x denote a geodesic from x to y of length at 
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most r, and j' x (0) its tangent vector at x. By construction, 7^(0) G C(w x ,[3 x ) and hence it forms 
an acute angle with v. This means that their images form an acute angle in the normal coordinates 
given by the exponential map at x. 

Now consider the normal coordinates given by the exponential map at y. In these coordinates 
7^ is a radial straight line emitted from the origin. Gauss's lemma (see [T7]) tells us that the angle 
between 7^ and j x in the normal coordinates at y is acute if and only if the angle between them in 
the normal coordinates based at x is acute. We have already shown that this second angle is acute. 

This means that in the normal coordinates given by the exponential map at y, j v must remain 
inside B(0, r) for some positive amount of time. Since this is true all x we know that the integral 
flow does not leave L r . Furthermore, the integral flow of X is always traveling towards K as it lies 
in W' . For each x S Kt r _g jr+ gi , let X x be the rate at which the integral flow of X at x is traveling 
towards K. Since X x > for all x £ K[ r _g >r+ g\ and K^ r _ s ^ r+ ^ is compact there is some A > which 
forms a lower bound on how fast the integral flow of X travels towards X. 

Construct the deformation retraction from L r to K r ^g by following each point along the flow 
of X until it reaches K r _g and then remaining stationary. The uniform lower bound on how fast 
the integral flow of X travels towards K combined with the observation that every point in L r is 
at most 25 from K r _g, tells us that in a finite amount of time every point in L r will be sent to one 
in K r _s- □ 

Theorem 4.4. Let fi G (0,1), r > 0. Let M be a smooth manifold with nowhere negative cur- 
vature such that every point has an injectivity radius greater than r. Let L a compact subset with 
dn(K,L) < 5. Suppose that there are no ^-critical points in K^ r _s iT -s+AS/fi] an d (4 + H 2 )o~ < ^ 2 r. 
Then L r deformation retracts to K r _g. 



Proof. Lemma 4.1 (with a = r — 5 and R = r) ensures that there exists a minimal cone field W over 



Kt r _g )r+ fii for Kg with length r. By Lemma 4.3 it is sufficient to show that there exists a Lipschitz 
vector field strictly subordinate to the complementary cone W of W . We will construct a useful 
cone field which is larger than W which has good continuity properties. For every cone field V such 
that V > W we know that W > V and hence it will be sufficient to find a Lipschitz vector field 
strictly subordinate to V . 

Since (4 + /i 2 )<5 < [i 2 r there is some ei > such that (4 + /i 2 )<5 < [i 2 r — 4ei, and hence R := r + ei 
satisfies 

5<{r-8-e 1 )- A T ^- 2 R. 

There also must exist some €2 > such that the injectivity radius is greater than R = r + £2- Since 
there are no /i-critical point in K[ r _g^ r _g +4: g/ C } and ||Va-(x)|| is a lower semicontinuous function of 
x, there must be some €3 > such that there are no /i-critical points in Kt r _g_ €3tr _g_ €3+ Mg +e3 \/Qj. 
Take e to be the minimum of these e^. This e will satisfy all three properties. 

Since M is a smooth manifold, it is triangulable. Barycentrically subdivide until we have a 
triangulation whose simplices have diameter at most e. Let £ be the simplicial complex by taking 
the closure of the union of all the simplices whose intersection with K\j._g^ r+ g^ is non-empty. From 
the bound on the diameter on the cells in £ we see that 

K[r-8,r+8] C S C K[ r _g_ e>r+ g +e y 

For each x G £ let 

Y(x):=l |J B(y,r)\nK 5 . 
\{y&A\xeA} I 
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From the restriction on the diameter on any cell in the triangulation and using the triangle inequality 
we know B(x, r) n K$ C Y{x) C B(x, R) n K$. Given any x G E there is some cell A containing 
x such that A n K\ r _^ rJr g\ is nonempty and hence contains some point y. This point y satisfies 
B(y, r) n Kg ^ which forces Y{x) to be non-empty. 

(with a = r — 5 — e and R = r + e) tells us that for every x G 2£V r _ < 5_ e ^a^i there 



Lemma 



4.1 



exists a unique acute minimal spanning cone of length R for K§ from x. Since Y(x) C A"^ n i?) 
and is not empty, we conclude that there is also a unique minimal spanning cone C(u x ,a x ) 

for Y(x) from x of length R. 

Set to be the cone field with attached cones C(u x , a x ) and V to be the complementary cone 
field of V. When we restrict V to K\ r _g r+ g\, our construction process ensures V > W and hence it 
is sufficient to find a Lipschitz vector field strictly subordinate to V . 

Given A C E, every x G A° has the same ^(x). Let xa be some point in A°. Observe that 
Y(x&) is compact. Construct V to be the minimal cone field over A for Y(x&) of length R. This 
minimal cone field exists because for any y G A we have B(y,R) n 3^(;ea) and ^(xa) C Y{y) 



(and as observed before there is a minimal cone from y for ^(y) of length R). Proposition 2.2 (with 
L = A and ii" = tells us that V is Lipschitz over A. Since V and V" agree in A° we know 

that V is Lipschitz over A°. This implies that V is also Lipschitz over A°. 

For every path 7 : [0, 1] — >• A with 7(1) G dA with j(t) G A° for t G [0, 1) we have that 
limt_s.i y(7(t)) C Y{ n f{l)) which implies that 

limC(n 7 ( t ),7r/2-a 7(t )) 5 C(it 7(1 ), vr/2 - a 7 (i)). (4) 

Informally we could say that that V' has some kind of lower semicontinuity. 

We will construct our unit vector field X over E which is strictly subordinate to V inductively 
over the skeleta of S of increasing dimension. First assign to each 0-simplex in £ any unit vector 
in the interior of the the cone in V' at that point. Note that this is always possible as the angle of 
the cone is always positive. 

Now suppose, inductively, that we there exists a Lipschitz unit vector field, Xf~, over the k - 
skeleton of E which is strictly subordinate to V . Let A be a simplex in E of dimension k + 1. By the 
inductive hypothesis we have some vector field X^q^ defined over dA which is strictly subordinate 
to V . We want to extend X^ to a Lipschitz vector field over A. 

Since V is Lipschitz on A° it has a unique continuous extension V" over A, which is also 
Lipschitz. From Q, the cone in V' at each point on <9A,is contained in the corresponding cone 
of V" . This implies that Xk\dA is strictly subordinate to V"\qa- Denote the cone field V" by 
{(x,C(u x ,a x ))}. 

Set E to be the fibre bundle over A where at each x we assign the subset C(u x ,a x )° of the 
unit tangent sphere. We can consider X^q^ as being a section of E\g&. The fibre bundle E has 
contractible fibers and a base space with trivial fundamental group so there is no obstruction to 
extending X^q^ to a continuous section of E defined over all of A. 

The space of Lipschitz vector fields which are X^ when restricted to dA is dense in the space 
of continuous vector fields which are X^ when restricted to dA. Furthermore, continuous sections 
of E which are X^ when restricted to dA form an open subset of continuous vector fields over A 
which are X^ when restricted to dA. This means there must be a Lipschitz vector field suitably 
close to Xk + i which extends X^. Set X^+i to be this Lipschitz vector field. 

If A is a k + 1 cell then X^+i is strictly subordinate to the corresponding V" . Since V and V" 
agree on A° we know that X^ + i is strictly subordinate to V . Thus X^i is strictly subordinate 
to V in the interior of every k + 1 cell. Since X^+i is also strictly subordinate to V' on the the 
/c-skeleton we can conclude that X^+i is strictly subordinate to V on the (k + l)-skeleton. 
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Inductively we can construct a Lipschitz vector field X over £ strictly subordinate to V . By 



our construction X is strictly subordinate to W' which allows us to apply Lemma 4.3 to build a 
defamation retraction from L r to K r _$. 

□ 



By doing the same process but using Lemma 4.2 instead of Lemma 4.1 we get the analogous 



theorem for when the ambient space has its sectional curvature bounded below by —1. 

Theorem 4.5. Let fi G (0, 1), r > 0. Let Ai be a smooth manifold whose sectional curvature 
bounded below by —1 and with an injectivity radius greater than r. Let K and L be compact subsets 
with dff(K, L) < 5. Suppose that there are no \i- critical points in K\ r _s,r-S+UI and 95 < 2tanh(r— 
5)fi 2 . Then L r deformations retracts to K r _g. 

5 Applications to point cloud data 

In this section we now consider the situation where we have some unknown compact set A which 
we are wanting to understand and we can sample A to generate a (potentially noisy) point cloud 
of A which we will denote by S. Historically conditions geometric conditions have been given on 
A for when offsets of the S and A are homotopic. This is because A is often assumed to have nice 
geometric properties whereas S, as a point cloud, has many critical points of its distance function 



nearby. The corresponding theorem produced using Theorem 4.4 and Theorem 4.5 with K = A 
and L = S is as follows. 

Corollary 5.1. Let ^ £ (0, 1), r > 0. Let Ai be a smooth manifold with sectional curvature bounded 
by k and whose injectivity radius is greater than r. Let A be a compact subset of Ai and S be a 
(potentially noisy) point cloud of A. Suppose that there are no ^-critical points in Ai a u. Then S r 
deformation retracts to A r whenever 

dn(S, A) < min /r — a, — A — \ , and dn{S, A) < ^ — ^ if k = 

or 



d H (S,A) < min <{ r - o, ^ — ^ } , and9d H {S,A) < 2tanh(r - d H (S, A))^ »/« = -!. 



4-/x 

Furthermore if A r _s deformation retracts to A then S r deformation retracts to A 



Proof. In order to apply Theorem |4.4| or Theorem 4.5 we need to make sure that 



[a, b] D [r-d H (S,A),r-dH(S,A) + 4d H (S,A)/ri 

and also that dn(S,A) < or 9du(S,A) < 2tanh(r — djj(S, A))fj? respectively. □ 

Of general interest is finding the homotopy type of A rather than A r . However, a sufficient 
condition for At, to deformation retract to A a is that there not be any 0-critical points in Ab\A a 
|15j . Also there are many shapes, such as like hairy objects, for which many offsets have a defamation 
retract even if there are small 0-critical values. 

We now want to present a paradigm for finding sufficient conditions on point cloud data for 
reconstructing any compact subset, lying in any Riemannian manifold, which has positive weak 



feature size. The first observation we need is that for Corollary 5.1 it is sufficient to have lower 
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bounds on the sectional curvature and the injectivity radius only for the points in A§ r and A^ r 
respectively. This is because no points outside this region are use in any of the proofs. Since A is 
compact there is some r > such that the injectivity radius of every point in A^ r is greater than r. 
Reduce r if necessary to ensure that r < wis(A) where wfs(A) is the weak feature size of A which 
we have assumed is positive. A% r is compact so there is some finite lower bound on the sectional 
curvature for points in A% r . By rescaling the metric on the ambient manifold if necessary (and with 
it scaling r) we can assume that the lower bound on sectional curvature is or —1. This means we 
can apply Corollary |5.1| It is clear a suitable /t and bound on da {A, S) in the Corollary must exist. 
Because r < wfs(^4) we can further state that the S r deformation retracts to A. This paradigm 
of reconstruction processes shows that what the ambient manifold is does not pose a theoretical 
barrier to the existence of reconstruction proofs. 

An alternative approach, as pointed out in [7J, is to consider geometric properties of S (or in 
their case offsets of S) itself rather than A. We can take his approach because S is a compact set 
and we do not require any smooth structure. This means we can also conclude another corollary 
with K = S and L = A. 

Corollary 5.2. Let [i £ (0, 1), r > 0. Let M be a smooth manifold with sectional curvature 
bounded by n whose injectivity radius is greater than r. Let A be a compact subset of M. and S be 
a (potentially noisy) point cloud of A. Suppose that there are no fi-critical points in S[ a ^. Then S r 
deformation retracts to A r whenever 

d H (S, A) < min jr - a, ^ — — \ , and d H (S, A) < ^ T 9 if k = 

or 



H J " v ' ' 4 + /x 2 



dn(S, A) < min Ir — a, — > , and dn(S, A) < - tanh(r — dn(S, A))fx 2 if k 



Furthermore if A r is homotopic to A then S r is homotopic to A. 

When the ambient space is Euclidean, it is reasonable to want to compare our reconstruction 
process to previous ones in the literature. Since these have been quantified in terms of //-reach we 
can first compare the required Hausdorff bounds on 5 := dn{A,S) where A is a compact set with 
^-reach > and S is a point cloud. If we consider the limiting case when r — 5 + 45/C < we 
can apply our reconstruction theorem (in the case of k = 0) once S/r^ < /i 2 /(4 + 4//). Notably this 
is an improvement on the bounds presented in [7], which is S/r^ < n 2 /(5fj, 2 + 12), for all \x and an 
improvement on the bounds in [2], where it is 

5 -3n + 3/i 2 - 3 + X/-8/U 2 + 4/x 3 + 18/i + 2/i 4 + 9 + ^ - 4/x 5 



r M 7fi 2 + 22/x + - 4/i 3 + 1 

for fj, < 0.945. 

One advantage of the approach of this paper is not having any requirements about the absence 
of //-critical points very close to A. A severe limitation of restricting to set with positive /z-reach 
is the inability to cope with sets that have cusps. At cusps the the /i-reach is zero for all values of 
ji > 0. The method used to overcome the shortfalls of /i-reach is to consider offsets of the compact 
set. For a compact set K, there are no //-critical points of dx in Kt a u if and only if the /x-reach 
of K a is at least b — a. This means we can compare different reconstruction theorems in terms of a 
lack of fi critical points in an annular region. 
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Let us assume that there are no //-critical points of dx in Kuu • For our reconstruction process 
we need 

J [i(b — a) fi 2 b 
\ 4 ' 4 + 4/t 



6 < min 



Here we would use r = . In comparison, for the reconstructions in [7] would need 

(6 - a)M 2 
5// 2 + 12 

and the the reconstructions in [2] we would need 



. , -3/i + 3// 2 - 3 + + 4/t 3 + 18/i + 2// 4 + 9 + // 6 - 4/z 5 
^ <(6 " a) 7m 2 + 22^ + ^-4^ + 1 

which is significantly worse when b — a is small in comparison to b. 

One possible future direction is the study of persistent homotopy. Persistent homology has been 
around for decades as a method of analyzing data [11] . It involves computing the homology classes 
of different offsets and studying when different homology class are born and die. A similar process 
could be done for homotopy. For a < b there is an inclusion map i : A a — > which induces maps on 
the homotopy groups t* : TTk{A a ) — > tt^A^) with changes only when a critical value of the distance 
function is passed. This paper gives sufficient conditions for S a > and Sy to deformation retract to 
A a and A^ respectively, and hence the map induced by inclusion nk(S a i) — > iTk(Sy) is the same as 
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